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The recent discovery of the putative 125-GeV Higgs boson has motivated a number of attempts 
to reconcile its relatively large mass with the predictions of the minimal supersymmetric standard 
model (MSSM). Some approaches invoked large trilinear supersymmetry-breaking terms At between 
stops and one of the elementary Higgs fields. We consider the possibility that electroweak symmetry 
breaking may be triggered by supersymmetry breaking with a large At, large enough to generate 
a composite field with the same quantum numbers as the Higgs boson and with a non-vanishing 
vacuum expectation value. In the resulting vacuum, the usual relation between the gauge couplings 
and the Higgs self-coupling does not apply, and there is no reason to expect the same upper bound 
on the mass of the lightest Higgs boson. In a simple model where the bound state is assumed to 
have no mixing with the other fields, we calculate the critical coupling At necessary for symmetry 
breaking using the lowest-order Bethe-Salpeter (BS) equation. Study of the BS equation is com- 
plicated by the structure of its lowest-order kernel, which is a crossed box graph, but we find an 
accurate approximation to its solution. In a realistic model, the mixing of the bound state with 
the fundamental Higgs boson creates a symmetry-breaking seesaw. We outline the steps toward a 
realistic model. 

The recently discovered 125 GeV boson [2 [2] is widely thought to be the Higgs boson. The mass outside the 
range predicted by the simplest supersymmetric models, combined with the lack of evidence for superpartners [S], 
has encouraged a number of efforts to reconcile low-energy supersymmetry with a relatively heavy Higgs boson (see, 
e.g., Refs. [^13). The models usually assume heavy masses for the superpartners, as well as some novel features, 
for example, strong couplings in the supersymmetry breaking sector [5]. While many models with gauge-mediated 
supersymmetry breaking predict a small tri-linear supersymmetry breaking coupling A, a large value of such coupling 
is, in fact, helpful in raising the range of the Higgs boson masses toward 125 GeV f4ji7j. Large tri-linear terms can 
appear in gauge-mediated supersymmetry breaking models, albeit some fine-tuning of parameters may be required in 
a realistic model jl]. 

However, if the tri-linear couplings are large, the low-energy realization of supersymmetry may differ dramatically 
from the usual set of predictions. It has been pointed out [HE] that the exchange of the (lighter) Higgs boson between 
(heavier) squarks can lead to formation of bound states, resonances, and a new strongly coupled realization of the 
minimal supersymmetric standard model (MSSM). Here we reconsider this possibility and will focus, in particular, on 
the possibility that supersymmetry breaking may trigger electroweak symmetry breaking via formation of relativistic 
squark bound states having the quantum gauge numbers of the Higgs boson. Such new states can mix with the 
Higgs boson, they can acquire a vacuum expectation value (VEV) , and the resulting multi- Higgs low-energy effective 
theory may have a very different appearance from the usual weakly coupled MSSM. At the same time, the ultraviolet 
behavior of the theory is preserved, and supersymmetry provides the usual solution to the hierarchy problem. The 
difference is in the low-energy effective theory, which contains different degrees of freedom: fewer squarks, and more 
Higgs bosons, whose VEVs produce a more complicated vacuum. In this vacuum, the usual MSSM relations between 
the gauge couplings and the scalar self-coupling do not hold, and, therefore, there is no reason for the upper bound 
on the lightest Higgs boson to be the same as in the usual version of MSSM. 

Let us consider a simplified version of MSSM, in which we will focus only on the third generation of squarks and 
will assume that only one tri-linear term is large: 

ASl-<l})iR + h.c. (1) 

where is the Y = 1/3 stop doublet under SU{2)l, if/ the Y = 4/3 stop singlet, and (j) the Y = -1 Higgs doublet. 
We omit writing cj)"^ terms. For simplicity we assume the squarks have a common mass M ^ a few TeV, considerably 
larger than the Higgs mass m. 

We have suppressed the SU(3) indices in Eq. [T] and will concentrate on the color-singlet bound state. This is 
the only bound state that can have a mixing with the fundamental Higgs boson, and as discussed below, we expect 
there to be a range of parameters in which this bound state has a non-zero VEV, while all the SU(3) non-singlet 
bound states (which can form through the Higgs exchange as well) have a zero VEV. This case corresponds to the 
standard-model-like multi-Higgs vacuum consistent with the data. 
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We seek a CP+ scalar doublet with Y = 1 (the quantum number of one of the Higgs fields) arising as a {ijit\^) 
bound state described by a Euclidean BS equation as shown in Fig. [l] 
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FIG. 1. The Bethe-Salpeter equation in "vertex" form. The double line represents the bound state $, of momentum P, and 
the single lines represent the constituent squarks, of mass M . 

This BS equation is in vertex form, where the internal lines represent propagators and the usual BS wave function 
is related to Q{P,p) by 



M 



^iP,p) = QiP,p). 



(2) 



It is useful to state the BS equation in this form because it is closely related to a gap equation whose non-trivial 
solution yields an estimate of how large the coupling in the kernel K must be to yield a bound-state Higgs that mimics 
the Higgs field </>. This bound-state Higgs with symmetry breaking mixes L and R stops and contributes to their mass 
difference. 

Let the line labeled {P/2) +p represent an outgoing f^, and the line labeled {P/2) —p represent an outgoing A 
few minutes of drawing Feynman graphs shows that the lowest-order kernel must be a crossed box graph, as shown 
in Fig.[2l 
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FIG. 2. The lowest-order kernel for the Higgs-like bound state. 



The lines labeled L,R are the stops, and the dashed lines are the Higgs fields of the MSSM. The Euclidean BS 
equation is 
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K{P,p,k). 



(3) 



We have omitted writing an SU (2) spinor index on Q and a corresponding factor Sij on the kernel. 

The complications of this one-loop kernel have prevented us from studying this BS equation at general momentum P, 
which would furnish a relation between the values of P for which the equation is solvable and the coupling At- Instead, 
we look for the value of At at which P = 0, corresponding to having four degenerate massless bound states. This 
total of four massless states corresponds to the two complex elements of the bound-state SU (2) spinor. Alternatively, 
these states can also be considered as states of a broken SU{2) x C/(l) theory: a zero-mass composite Higgs boson, 
plus three zero-mass Nambu-Goldstone bosons. We are interested in Higgs bosons whose mass is comparable to m, 
the lightest Higgs boson mass in the MSSM. By hypothesis this is close to 125 GeV, while the squark mass M is much 
larger. Thus, it should be a reasonable approximation to consider the bound-state Higgs as having zero mass and 
study the P = BS equation. In addition, this P = BS equation can be tackled with decent quantitative accuracy. 
From now on we use the notation K{p,k) = K{k,p) for the original kernel at P = 0, whose bound states can only 
occur for specific values of At that are eigenvalues of the BS equation. 

In the general case with P ^ the BS wave function ^'(P, p) is the Fourier transform of the matrix element 



(o|r(tli 



Xl)iRix2))\P) 



(4) 
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FIG. 3. Diagrams for the effective potential quadratic in SM^ (weights not shown). 

with P conjugate to the center-of-mass coordinate X = {l/2){xi + X2) and p conjugate to xi — X2- the state \P) 
is the bound state. At P = this looks hke a vacuum-to- vacuum propagator Ai/j(a;i — ^2), but one must be 
careful about what the vacuum means. Just as in superconductivity the true vacuum is a non-perturbative construct 
quite different from the bare vacuum; in our case the true vacuum has matrix elements connecting L and R stops. 
This connection comes from a symmetry-breaking order parameter that is a mass splitting SM'^{p) found in this LR 
propagator, vanishing in the symmetric case, that mixes L and R squarks. This is analogous to the ('0V') propagator 
of superconductivity |10j . To lowest order in this order parameter the diagonal propagators of the t^ jj fields are just 
those already shown in the BS equation: 

ALL{p)^AnBip)^-^^ (5) 

while the LR mixing propagator is 

Away from the symmetry-breaking threshold, one would expect that for larger values of At the mass of the Higgs 
particle moves away from zero, while the Nambu-Goldstone fields remain massless. This can proceed (for example, 
see [TT]) through a tachyonic solution to the BS equation, much as in the stabilization of a Mexican-hat potential 
where the stable non-perturbative vacuum yields a condensate and a Higgs boson of normal mass. As for the Nambu- 
Goldstone bosons, since the pioneering work of Nambu [101 we know that these massless Nambu-Goldstone excitations 
occur as a consequence of a non-trivial solution to a gap equation, an integral equation whose solution is a symmetry- 
breaking order parameter. This gap equation is essentially the BS equation at P = 0. (See [12] for a proof of 
this Nambu-Goldstone mechanism in gauge theories such as QCD.) We express the dynamics of symmetry breaking 
through the usual two-particle irreducible (2PI) effective potential F T3], in which F is a functional of (SM^(p^). In 
this first investigation of the bound-state Higgs we ignore a number of interesting phenomena, including the VEV of 
the elementary Higgs fields and their possible mixing with the bound-state Higgs, so the effective potential (in the 
notation of [13]) is 

F = ^TrjlnG + [1 - GG^^]} + 2PI graphs, (7) 

where the trace is over space-time as well as other relevant indices, such as particle type, G is the exact propagator, 
and Go is the free propagator (when relevant; the term in square brackets is omitted for the LR propagator). The 
extrema of F as the G are varied yield the Schwinger-Dyson equations of the theory. 

To lowest order in SM'^ the effective action is given by the diagrams shown in Fig. [sj which give for F the expression 

^ d^pp{p^)[5M'{p^)]^ - f d^p f d^kp{p')SM\p^)K{p,k)p{k^)SM\k') (8) 



2 J '^'^'■^ " "-"^ " 2(27r 

where 



Variation of the quadratic terms in F with respect to (5M^ yields 



mp)' = 7^ f d'k-^sM{kr-i—Kip,k). (10) 
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This equation is analogous to standard gap equations for chiral symmetry breaking. Just as for chiral gap equations 
it is in fact the original BS equation, in vertex form, at P = 0, illustrating as before [101 [12] the necessary existence 
of composite Nambu-Goldstone bosons when symmetries are broken without elementary Higgs fields. It differs from 
chiral symmetry breaking gap equations because the kernel is well-behaved in the UV and there are no UV divergen ces . 
The kernel falls like (modulo logarithms) at large momentum, implying the same falloff for JM^, and Eq. (10) 
is finite. 



To analyze Eq. (10 1 we need to analyze first the kernel K(p,k). This kernel has the form 

1 



1 



(27r)4 
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+ A'P] [{p + lY + m2] [{k + l-pY + Af2] [{k + If 
1 

Ifi 



(11) 



with 



D = k'^{xiX2 + ^3X4) +p^(xia;4 + X2X^) + (p + kYx2Xi + {p - k fxix^ + M'^{x2 + .T4) + m'^{xi + X3). 



(12) 



Now suppose that Af^ ^ m^, in which case X2,X4 have to be small compared to the other Feynman parameters. 
So write X2 = Ax, 2:4 = A(l — x), with new integration variables running from to 1. The integral over A will be 
dominated by small A, so we can drop this variable judiciously. Then approximately 



and 



where 



dxiS ^1 — Xi^ = XdXdxdxidx3S{l — xi — X3) 

D = xixsip - kf + rrf{xi + Xj) + \[ak'^ + (1 - a)p^ + Af^] 
a = XiX + 2:3(1 — x), < a < 1, 



and we dropped a term ^ A^ in D. 

Now do the integral over A explicitly, with the result 



Kip,k) 



1 

167r2 



dxidx3dxS{l — xi — X3) 



In 



A + B 
B 



1 

A{A + B) 



where 



A = [ak^ + (1 - a)p^ + M\ B = XiXz{p - kf + nf. 



(13) 
(14) 
(15) 

(16) 
(17) 



We now show that at large k^ ^ p^ , K ^ (1/fc^) Infc^ and so vanishes rapidly. This means that the integral in the 
gap equation also vanishes rapidly, as we will soon see. Since B does not depend on x the integral over x can be done, 
yielding 



K 



lnfc2 

87r2/c4' 



(18) 



a result that we can get exactly at p = 0. 

The next step is to reduce the gap equation to a one-dimensional equation by integrating over the angles oi k. In 
the gap equation, because the angle between the four-momenta appears only in logarithms or in a term parametrically 
small with respect to Af^, we make the approximation 

j dVLkF[{p ~ kf] « 2^\e{p' - e)F{p') + e{k^ - /)F(fc')]. (19) 

This is exactly true for F — \/{p~kf or for constant F, and is acceptable for the logarithmic functions we encounter. 
After projecting out the s-wave, to lowest (quadratic) order in SM'^ the relevant part of F is 

„2 „2 „/ 2Nrr7i,r2/ 2m2 



dp'p'p{p')[5M\p')Y 



2(27r)4 



dp'p' / dk' e p{p')5M'{p')K{p\e)p{e)5M\e) 



(20) 
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where K is the s-wave projection of the kernel. 

Variation of this equation yields the s-wave projection of Eq. (10), which becomes a standard one-dimensional 
homogeneous Fredholm integral equation with a discrete spectrum of eigenvalues . We seek the lowest eigenvalue 
by inserting a trial function into Eq. ( 20 ) and doing the integrals numerically, including an approximation to the 
integral over xi in the kernel (see Eq. (13)), for various values of the Higgs-to-stop mass ratio m/M. 

We motivate our trial functions for the crossed-graph kernel of interest from known exact results [14] for the BS 
equation with the massless kernel 



K 



(p - fc)2 • 

For the vertex form of the BS equation, given in Eq. ([3|, the lowest eigenfunction at P — is 

1 



QiO,p)=dM\p') 



p2 + ' 



(21) 



(22) 



Naturally, the class of trial functions of this form with AP replaced by a variational parameter yields the exact 
result. In the present problem the asymptotic behavior is different, so we choose as a zeroth-order trial function 



1 



(23) 



We have studied other trial functions, such as l/(p^ -I- /i^)^, with similar results. We improve this first variational 
estimate by using SMq as input to the right-hand side of Eq. (20), numerically calculating a new output SMi- We 



made a simple but accurate fit to SAd^, amounting to adding a term ~ p'^ to the denominator of Eq. (23). Then we 
used the average (5M| = (1/2)((5Mq -\-5M1) as a trial function, and calculated the output again. This yields excellent 
agreement between the new input and output, as shown in Fig. |4j for the specific value m/M — 0.05. 




FIG. 4. A comparison of the input and output using the second-order eigenfunction 5M2 , as a function of p /M 
numerically as described in the text. In the case of the exact solution, the two curves would be identical. 



calculated 



At this elementary Higgs mass the critical coupling resulting from our numerical calculations is 

15.14. 



M 



(24) 



This estimate is to be compared with the value of At/M = a/6 that maximizes the lightest Higgs mass in an 
approximate one- loop calculation [4|[T5]. But this is not the final verdict, since mixing of the bound Higgs with the 
MSSM Higgs and other possible bound states need investigation. 

Also of interest is the needed critical coupling for various masses of the elementary Higgs field. This is shown in 
Fig. [5] As expected, the critical coupling increases with increasing Higgs mass. 
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FIG. 5. Behavior of the critical coupling as a function of the elementary Higgs mass. 



While it appears plausible that supersymmetry breaking in the MSSM can trigger electroweak symmetry breaking 
via the formation of bound states with non-zero VEVs, further work is needed before one can build a realistic model 
and compare its predictions with the data. In addition to the color-singlet states, the same trilinear scalar interactions 
can cause colored bound states to form. The viability of the model depends on its ability to produce a standard-model- 
like vacuum with broken SU(2)xU(l) but unbroken SU(3), in which case the colored bound states are turned into 
color singlets by strings attached to gluons or quarks. It is well known that the MSSM, in its traditional realization, 
has a number of dangerous color and charge breaking minima, although cosmological evolution favors the vacuum with 
unbroken SU(3) even in some cases where it is not the global minimum of the potential [TB]. In our case, one must now 
re-examine the same issue taking into account a number of new effective degrees of freedom. While the full analysis 
is obviously complicated and the results will inevitably be model-dependent, there is one feature of the color-singlet 
states that sets them apart from the rest. The color-singlet states can have a mixing with the fundamental Higgs 
bosons via the same coupling as that which enters the BS equation. The mass matrix in the bound-state-Higgs basis 
has both diagonal terms and off-diagonal terms. In contrast, the colored bound states can only have diagonal terms. 
As the BS coupling increases, the mass squared of each bound state decreases, as discussed above. Since the scalar 
exchange forces are essentially color-blind, the bound states with different SU(3) properties can have similar binding 
energies. However, thanks to the off-diagonal terms, the colorless bound states can develop a VEV simultaneously 
with the Higgs boson for some value of the trilinear coupling for which the diagonal terms are still positive. This 
possibility leads to an appropriate standard vacuum. 

The mixing of the fundamental Higgs and the bound state occurs through a diagram involving the solution to 
the BS equation, as shown in Fig. [6] Also shown is an approximation to the four-point coupling of the composite 
Higgs, again expressed in terms of the BS solution. This quartic coupling is not related to the gauge couplings by 
supersymmetry because it depends explicitly on the supersymmetry breaking parameter At . Therefore, in the vacuum 
where the fundamental Higgs boson and the bound state mix, the usual upper bound on the lightest Higgs boson 
mass does not apply. 

Electroweak symmetry breaking occurs via a symmetry-breaking seesaw^ when the following mass matrix has a 
negative eigenvalue: 



Mt 



(25) 



Here /i^ is the Higgs mass-squared parameter, which is negative in the Standard Model, but which we take to be 
positive; A/hb is the mixing parameter calculated from the diagram shown in Fig.|6j and Afe is the mass of the bound 
state. Our calculations above dealt with M| = = 0, which required the large value of At/M quoted in equation 



(24). However, in the presence of mixing, the symmetry breaking can occur for a non-zero Mb and a positive /i^, as 



long as 



detTWn = Mh^-^b 



A^B 



< 0. 



(26) 





FIG. 6. Mixing of the composite and elementary Higgs fields, and four-point couplings of the composite Higgs field. 



The value of the mixing parameter Mhb can be extracted from a solution of the BS equation with a proper wave 
function normalization. We leave this calculation for future work. Here we parameterize 



HB 



(27) 



where ^ is dimensionless number (which may depend on At and on other parameters). Now the critical value of At is 
different from that in Eq. [24] 



M 



15.4 



(28) 



Given that can be arbitrarily close to zero (and, in fact, /i^ < in the Standard Model), the symmetry breaking 
can occur even for smaller values of At/M . Of course, if /.t^ has to be fine-tuned to be small, the scenario in question 
becomes rather contrived. For reasonable values of parameters and for the binding energy that is not as large as the 
mass, which is the case for a smaller A^, there is a self-consistent set of values, for example. 



- 10^ GeV, Mb 2M - 2 TeV, At 



(i^Mb 300 GeV 



and, thus, 



A 

M 



0.3 



(29) 



(30) 



One can now ask whether these values are consistent with vacuum stability with respect to dangerous color and 
charge breaking minima in the potential [ini [T7]. In general, one has to analyze the full potential of the MSSM to 
identify the global and the local minima. If the global minimum is color and charge preserving, the model is viable. 
If the standard-model-like vacuum is only a local minimum, one has to examine cosmological evolution because in 
many cases the vacuum that is populated first is the one in which the squarks have zero VEV [TBI HH] • One then 
has to calculate the lifetime of the false vacuum and determine whether the universe is likely to undergo cooling and 
remain in the false vacuum until the present time. An approximate criterion for the absence of dangerous color and 
charge breaking minima, based on the analysis of Ref . |16j , is 



A 
M 



^ 2.7, 



(31) 



which can be consistent with Eq. (30) in some reasonable range of parameters. 



A large trilinear coupling at low energy can become even larger at higher scales, but the effects of renormalization 
group running depend on the size of the coupling. For negative sign, the initially large At can cross zero at a high 
energy scale |4]. 

In addition to the SU(2) doublet bound state, the large At can lead to formation of a gauge singlet state in which 
{II, t*j^} in the initial state exchange the Higgs boson and convert into {tji, which, in turn exchange the Higgs 
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boson and convert back to {II, t^}- This is similar to the doublet bound state we have considered, but with a 
BS wavefunction appropriate for the one-particle kernel. The existence of such a gauge singlet scalar is a generic 
prediction of our modelF] 

So far, we have focused on the version of MSSM in which only the top trilinear term is large. However, if the 
other trilinear couplings and the fj,B term are large as well, new additional states can form through the Higgs boson 
exchange, and new diagrams can contribute to the bound states we discussed. 

The electroweak precision measurements should impose constraints on any strongly coupled model built along the 
lines we discussed. These constraints imply a lower bound on the mass of the squarks. It is possible, and, in fact, likely, 
that the squark masses would have to be of the order of 5 — 10 TeV for the model to be consistent with the precision 
measurements. This forces the A term to be correspondingly larger, and there may appear to be a small hierarchy 
between the supersymmetry breaking scale and the electroweak scale. This hierarchy may impose some degree of 
fine-tuning on a realistic model based on strongly coupled broken supersymmetry. This is a likely potential drawback 
of the otherwise very appealing scenario, in which the scale of the electroweak symmetry breaking is determined by 
the breaking of supersymmetry. However, the class of models we discuss still possesses a robust solution to the big 
hierarchy problem: above the scale of the bound states, the MSSM exists in its usual incarnation, and supersymmetry 
stabilizes the scales in the usual way. 

We have examined the possibility of electroweak symmetry breaking by the formation of bound states of squarks 
via the Higgs boson exchange, which bound states can mix with the fundamental Higgs bosons and can acquire VEVs 
simultaneously with these fundamental bosons. This scenario is clearly different from the widely discussed technicolor 
models [TH] , including walking technicolor [50] , models with color singlets [2T] , and the models in which supersymmetry 
and technicolor are combined |22j . Our scenario has the potential to relate the scales of supersymmetry breaking 
and electroweak breaking in a new way, but the applications to realistic models requires a more detailed analysis. 
It is evident that bound states can play an important role, and they should not be neglected in models with large 
trilinear supersymmetry breaking terms. The next step that we plan to undertake is to investigate the mixing of 
the elementary and composite Higgs bosons and the resulting symmetry-breaking patterns. These results will be 
presented elsewhere. 
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